Dynamics of Relativistic Solitons 



Daniela Farinaf and Sergei V. Bulanov|§ 

f Istituto di Fisica del Plasma "Piero Caldirola" , Consiglio Nazionalc dcllc 
Ricerche, EURATOM-ENEA-CNR Association, Milan, Italy 
\ Advanced Photon Research Center, Kansai Research Establishment, JAERI, 
Kizu, Kyoto 619-0215, Japan 

§A. M. Prokhorov Institute of General Physics, Russ. Acad. Sci., Moscow, 
Russia 



Abstract. Relativistic solitons are self-trapped, finite size, electromagnetic 
waves of relativistic intensity that propagate without diffraction spreading. They 
have been predicted theoretically within the relativistic fluid approximation, and 
have been observed in multi-dimensional particle in cell simulations of laser pulse 
interaction with the plasma. Solitons were observed in the laser irradiated plasmas 
with the proton imaging technique as well. This paper reviews many theoretical 
results on relativistic solitons in electron— ion plasmas. 



1. Introduction 

Relativistic solitons are self-trapped, finite size, electromagnetic waves of relativistic 
intensity (eE/mujc > 1) that propagate without diffraction spreading. Many different 
physical effects play a role in the formation of relativistic solitons: dispersion effects 
due to the finite particle inertia, nonlinearities due to relativistic mass increase, as 
well as ponderomotive effects which force the plasma density redistribution. 

The theoretical investigation of relativistic solitons in electron-ion plasmas is a 
relatively old problem in plasma physics, which has been treated by many authors 
in the past, and has recently gained new attention in the literature ^ EI El El El El 
□ ElElEUOIllHIISmDjailEl- The analyses have been performed mainly in the 
framework of the ID relativistic fluid approximation, in which solitons are described 
by the solutions of a set of coupled nonlinear differential equations for the electrostatic 
and electromagnetic potentials with suitable boundary conditions. 

Relativistic solitons have been seen in multi-dimensional particle in cell (PIC) 
as well as in fluid simulations of laser pulse interaction with the plasma ^| ^1 1191 
1201 1211 17^1 I2HI I24| . These solitons consist of electron density depressions and intense 
electromagnetic field concentrations with a larger amplitude and a lower frequency 
than those of the laser pulse. Since a significant portion of the overall electromagnetic 
energy is trapped in the form of solitons, solitary waves can play an important role in 
the laser-plasma interaction. 



2. ID relativistic solitons 



Here, the theory of ID circularly polarized solitons is presented within the relativistic 
hydrodynamic approximation used to describe both the electron and ion components 
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|51 llUllTT] . The plasma is assumed to be cold with zero ion and electron temperature. 
Length, time, velocity, momentum, vector and scalar potential, and density are 
normalized over c/w pe , w pe , c, m Q c, m e c 2 /e, and n , respectively, being uj pe — 
(47rnoe 2 /m e ) 1 / 2 the electron plasma frequency, m a the rest mass with a = e, i, and 
no the unperturbed electron (and ion) density. In the Coulomb gauge, the Maxwell's 
equations for the vector and scalar potentials, A and <p, and the hydrodynamic 
equations (continuity, and momentum balance) for the densities n a and the canonical 
momentum P Q of electrons and ions can be written as 

d 2 d 

AA - —A - = n e v e - n i v l , (1) 

= n e - Ui , (2) 

d 

— n a + V • (n Q v Q ) = 0, (3) 
at 

d 

— P Q = -V(p Q + 7a ) +v a x V x P a , (4) 
at 

where P Q , and j a are related to the kinetic momentum p Q by P Q = p a + p a A, and 
7 Q = (1 + IPqJ 2 ) 1 ^ 2 with the parameter p a = (q a /q e )(m a /m e ), and v Q =p a /"f a is the 
fluid velocity. 

For the ID case in which d y = d z = 0, the vector potential is assumed of the form 
(circular polarized wave) 

A y + iA z = a(£) exp exp (— iwt + ikx) (5) 

with £ = x — Vt, while all the other quantities, 4>, n a , j a , and p xa , are assumed to 
depend only on the variable £. In the present case, the relations A x = and P^ = 
are satisfied, where the symbol _L refers to the direction of the laser pulse propagation. 

Imposing as boundary condition at the point £ = — oo, a = ±ao, (f> = 0, 
n a — 1, and p xa = (plasma at rest), the longitudinal component of the kinetic 
momentum, the energy, and the density of each species can be expressed as a function 
of the potentials as p xa = (Vip a — R a )/(l — V 2 ), ^ a = {ip a — VR a )/(l - V 2 ), 
n a = V{ j> a /R a - V)/(l - V 2 ), w here ^ a = T 0a + Pa <j>, R a = [V4 - (1 - ^ 2 )r 2 ] 1/2 , 
T a = yl-j- p 2 a 2 and To a — yl + p 2 of- Then, for an electron-ion plasma the 
following closed system of equations for the potentials is obtained 

d 2 4> v fi> e ^ , 

(6) 



d£? l-V 2 \R e Ri 

d 2(1 , ( -2 72 a o\ V ( 1 , p 



£ = -Mi-3) (8) 



aie ' " " a 4 J "'l-V 2 \R e ' R 

al6 

dt; \ a 

where Q = (w - kV)/(l - V 2 ), k = (k - wV)/(l - V^ 2 ), and p = = m^rn,. 

The system of equations (|6I7H describes coupled Langmuir and circularly polarized 
transverse electromagnetic waves. Note that the r.h.s. of (JJJ is proportional to 
n e he — Hi/7» and thus represents the balance between the striction nonlinearity due 
to perturbation of the density and the relativistic nonlinearity due to mass. Equation 
||SJ| describes the related phase evolution of the e. m. amplitude. For ao = 0, the 
solution of (JSJ is simply 8 = 9q — fc£, so that the e. m. field defined in JSJ has the 
simple form a(£) exp[— iujt + iOq], being r = t — Vx. For ao 7^ 0, the phase evolution 
is non trivial. 
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The system (|6I7|) can be put in Hamiltonian form, and has a first integral 

H (a, a', 0, 0') = i-^ fa' 2 + Cu 2 a 2 + #4 V^' 2 "^"- = *(9) 
2 \ a z ) I p 

where the symbol / denotes derivative with respect to £. The value of the constant K 
is determined by the boundary condition. 

For the case of a pure transverse electromagnetic wave (i.e., <fi = and a = an), 
from (JZJ) the dispersion equation is obtained 

u, 2 = n 2 + k 2 , n 2 = - 1 + - 9 (io) 

Q being the plasma frequency modified by the relativistic effects. Equation (|10|l 
corresponds to the Akhiezer-Polovin result |25) with the ion motion taken into account. 

So far, explicit reference to an electron-ion plasma has been made. However, 
putting p = 1, Eqs. Ij6l9|l are valid also in a electron-positron plasma |2fi| . 

The particular class of localized solutions of equations I|tjl7fl is investigated, 
satisfying a' = 0, a" = 0, <fr' = 0, and 4>" — at £ = ±oo. For sake of simplicity, the 
case k — is considered, so that k — ujV, and Co — u (see e.g., |S]). Note that here V 
is the group velocity. 

2.1. Quasineutral approximation 

Within the quasineutral approximation the main features of the localized solutions 
can be identified in a simple manner, assuming n e — ni, i.e., putting <f>" = 5 . Then, 
the electrostatic potential can be written in terms of the vector potential amplitude, 
<p = (ToiT e — ro e ri)/r, being r = Ti + pT ei and the equation for the vector potential 
in the quasineutral approximation becomes 

The above equation has the following first integral Ho(a,a') — 1/2 (a' 2 + <I> 2 a 2 ) + 
V v/r 2 - (1 - V 2 )T 2 /p(l - V" 2 ) 2 = K , and the solution for a is reduced to 
quadrature. Then, in the limit p 1, the following simple expressions hold: 
cj) w Vl + a 2 - V 1 + a 0' = rij w 1 + p(/)/V 2 , and w xe ^ w p^/V". 

From the analysis of the topology of the phase space (a, a') of the Hamiltonian 
H Q , different kind of solitary solutions are found varying the propagation speed, and 
are summarized in tabled 



Table 1. Solitons in a cold plasma. 



Propagation 


Solution 


Potential a 


Potential <j> 


Density 


speed 


kind 


at — oo, +oo 


profile 


profile 


V 2 > V s \ 


bright 


0, 


positive 


accumulation 


V 2 = v s \ 


shock wave 


0, oi 


positive 


compression 


V 2 < V 2 < V 2 


dark (gray) 


ao, ao 


negative 


evacuation 


V 2 = V 2 


shock wave 


a , 


negative 


rarefaction 


< V 2 < V 2 


dark (black) 


ao, -a 


negative 


evacuation 



Solitary solutions are found for oj 2 (1 — V 2 ) < 1 + p, with dark solitons satisfying 
the dispersion relation (|f Oj l . In both cases, the frequency spectrum corresponds to the 
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evanescent region in the (u>, k) space, i.e., to the region lo 2 — k 2 < 1 + p. When the 
ion dynamics is taken into account, bright solitons occurs for propagation speed larger 
than a threshold value scaling with */jo, and dark solitons occur below this threshold 
value. Standing bright solitons are found only when the ion dynamics is neglected. 
Solutions in the form of shock- waves are found at critical propagation speeds. These 
shock waves (occurring in the cold plasma limit) are collisionless, and similar to 
those found at ion-acoustic speed The characteristic velocities introduced in 

table n depend on the potential amplitude and have the following expressions in 
the low amplitude limit and for p -c 1: « p(l + of), V 2 « p(l + Oq/16), 

V 2 s» p(1 + Oq/2). In the same limit, bright solitons are described by a — a TO sech(K£), 
with k 2 = 1 + p — cj 2 (1 — V 2 ), and a 2 n = An 2 / '(1 — p/V 2 ), while dark solitons by 
a = aotanhOO, with a 2 , = An 2 /{p/V 2 - 1), and w 2 (l - V 2 ) = 1 + p- ag/2 01m. 

At large amplitude, quasineutrality is violated 0, and the set of coupled equations 
<|6I7|) has to be solved numerically. However, the above investigation on the existence 
of the different kind of solutions has a general validity. 

2.2. Bright solitons in e-i plasma 

In the case of large amplitude bright solitons, a wide class of solutions can be found. 
Here, the analysis is restricted to the class of solutions with single humped <f> profiles, 
and a profiles characterized by the number of nodes p. For fixed V, the eigenspectrum 
of the bright solitons can been determined solving numerically Eqs. I|6I7(I with suitable 
boundary conditions [SJ HOj . 

The eigenspectrum relevant for low p values is shown in figure ^ both for the case 
in which the ions are fixed, i.e., at p — 0, and for the case in which the ion dynamics 
is taken into account. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

V V 



Figure 1. Eigenspectrum of bright solitons as a function of the group velocity 
V. Cases (a) and (b) refer to p = (no ion dynamics) and p = 1/1836 (with 
ion dynamics), respectively. Both the regions corresponding to the continuos 
spectrum of single— peaked solitons and the discrete spectrum for p = 0, . . . , 4 are 
shown. The dashed curve in frame (b) delimits a region of no solutions. 

When the ion dynamics is neglected (i.e., p = 0), the bright soliton eigenspectrum 
has the following features. Single-peaked solitons (p — 0) admit a continuous spectrum 
for any velocity value In particular, standing stable solitons (V — 0) are found 
with frequency in the range 2/3 < lo 2 < 1, and max a < ^/3, as it has been shown 
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analytically in HJ. Multi-peaked solitons (p > 1) have a discrete spectrum for 
velocities larger than a small critical value, and at fixed V, the frequency ui decreases 
with increasing p. The peak value of the potentials a and <f> is found at the minimum 
velocity, and increases with p 

Quite different results are found when the ion dynamics is taken into account 
(i.e., p / 0), even at a finite propagation speed. Note that at high frequencies (and 
velocities), the ion dynamics does not play any role. Single-humped solutions are 
found for velocities larger than V s ~ ^fp 0, with both a continuous spectrum [T3] 
and a discrete spectrum ^Uj. For p > 1, only a discrete spectrum is found JU|- At 
low p values, the frequency is non monotonous as a function of V, and at V = Vbr the 
branch ends since the soliton breaks. The structure of the solution for p = 0,1,2 and 
a velocity value close to breaking is shown in figure [SJ where the potential waveforms 
are plotted. It is found that the ions pile at the soliton center, while the electrons 
at its edges, giving rise to very peaked density distributions, and breaking occurs 
because the ion density diverges in the center. At the same time, the ion and electron 
velocities tend to Vb r at the center of the soliton and at its edges, respectively, so 
that in these regions the particles move with almost the group velocity of the soliton. 
The ion velocity profile shows a cusp in the center, which is the signature of the 
nonlinear wavebreaking. The soliton breaking provides a novel mechanism for the ion 
acceleration in the high intensity laser pulse interaction with plasmas |1U| . 

Solitons at large p number have been investigated numerically and analytically in 
different approximations |rfl 171 1%1 IT3|. 
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Figure 2. Potential waveforms for bright solitons with p = 0, 1, 2 and velocities 
close to breaking. 



2.3. Dark solitons in e-i plasma 

As already pointed out, dark solitons are found only when the ion dynamics is taken 
into account and at very low velocities (V < \fp)- The solutions of the full nonlinear 
system for ao ^ are obtained numerically, as in the case of bright solitons • The 
obtained spectrum and the corresponding ao value are shown in figure |2| as a function 
of the velocity. At fixed V, the spectrum is continuous at low ao, and it becomes 
discrete above a threshold value, which increases with V. The discrete spectrum 
is made up by a large number of separate branches, each of them characterized by 
different densities and velocities profiles. The ion density profile show always a dip, 
while at large ao the electron density tends to peak in the region inside the soliton. At 
a critical amplitude value (dependent on V) breaking of the solutions occurs, due to 
the peaking of the electron density. Solitary solutions are found up to V rs 0.051, and 
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a 5.8 (for p = 1/1836). Above these values, no solutions have been found for the 
chosen parameters. Both black and gray solitons are found for V smaller and larger 
than V 8 , while a shock wave is found at V = V s . Figure 01 shows the waveforms of the 
potentials for dark solitons (black and gray) and the shock wave. 
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Figure 3. Eigcnspectrum of dark solitons (a) , and asymptotic c. m. amplitude 
ao (b) as a function of the group velocity V. The symbols GS and BS indicate the 
regions corresponding to the continuous spectrum of gray and black dark solitons, 
respectively. The curve corresponding to shock waves is also plotted. 




2-4- Warm plasma effects 

The results described in the previous section have been derived in the cold plasma 
approximation. However, it can be easily argued that solitons moving at "low" velocity 
can be affected by pressure terms. The theoretical investigation of temperature effects 
has been performed within the quasineutral approximation for standing solitons in 
a fully relativistic treatment |14) . and for the case of moving solitons in the weakly 
relativistic approximation JS]. The results are summarized in tabled 15 , which refer 
to small finite amplitude solitons and non relativistic temperatures. 

As expected, temperature effects play a crucial role in the range < V 2 < p + c 2 , 
being c s the ion-acoustic speed. Note that both inertia and temperature terms drive 
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Table 2. Solitons properties in a warm plasma. The following symbols are used: 
^te i = T ei i/m e ^c 2 , c 2 = Y e pV 2 e + r^V^, with V ej j the adiabaticity parameter. 



Propagation 


Soliton 


Potential <f> 


Density 


Species 


speed 


variety 


at the center 


at the center 


involved 


< V 2 < v?i 


bright 


positive 


evacuation 


i+e 


V?, <V 2 <c 2 


bright 


negative 


evacuation 


i+e 


c 2 s <V 2 <p + c 2 s 


dark 


negative 


evacuation 


i+e 


p + c 2 s < V 2 


bright 


positive 


accumulation 


i+e 


P + cl^V 2 


bright 


positive 


evacuation 


only e 



the process of soliton formation in a warm plasma. In particular, the formation of 
single-peaked solitons is observed at zero and very low propagation speed, lower than 
the ion acoustic speed. The relevant electron and ion density profiles show a dip in 
the center of the soliton. In the low temperature limit, density depletion in the center 
of standing solitons may become almost total |14|. 

3. Solitons in laser plasma simulations 

Relativistic solitons has been investigated with the help of fluid simulations and 2D, 3D 
PIC simulations for laser-plasma interaction [1711181 rO?ll21ll2*2ll23| . The mechanism of 
soliton generation by a high intensity laser pulse propagating in a underdense plasma 
can be descibed as follows. The laser pulse interacting with the plasma loses its 
energy, as it generates wake fields behind itself. Since the process is adiabatic the ratio 
between the e. m. energy density and the frequency is conserved, so that frequency 
downshift occurs for the wake field. Moreover, since the group velocity of the laser 
pulse decreases with the frequency, the low frequency part of the e. m. radiation of 
the pulse propagates with very low velocity. When the frequency becomes lower than 
the Langmuir frequency, the e. m. energy becomes trapped inside the related density 
cavity, thus forming a slowly moving solitary structure. Note that in 2D and 3D it is 
found that solitary waves can merge, thus they are not strictly speaking solitons. It 
has been observed that a large part of the laser pulse energy can be transformed into 
solitons |19j . It can then be argued that solitary waves may play an important role in 
the laser-plasma interaction. 

When PIC simulations are performed in inhomogeneous plasmas |2(J| . it is found 
that low frequency solitary waves generated by superintense laser pulses are accelerated 
along the density gradient towards the lowest density values (e.g., the plasma-vacuum 
interface) where they radiate their energy in the form of low frequency electromagnetic 
bursts. This process can be a signature of soliton formation. 

PIC simulations show that the time of the soliton formation is much shorter than 
the ion response time, so that ions can be assumed at rest during this process. For 
approximately (m-i/me) 1 / 2 oscillation periods of the e. m. field inside the soliton 
(which is of the order of 2ir/uj pe ), the ions can be assumed as fixed. For longer times, 
the ponderomotive force starts to dig a hole in the ion density, and the parameters of 
the solitons change. What was a soliton on the electron timescale is no more a soliton 
on the ion timescale. Nevertheless, a low frequency e. m. wave packet remains well 
confined inside a slowly expanding plasma cavity. This e. m. entity has been called a 
postsoliton [2*T] . 
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In addition, ID fluid simulations in "slightly" overdense plasmas (i.e., for 
plasma densities < 1.5 times the critical density) have shown that penetration of 
relativistically intense laser radiation occurs by soliton-like structures moving into 
the plasma 24 . 

4. Final remarks 

Many theoretical problems are still open even for ID solitons, as, e.g., the temperature 
effects at any propagation speed and arbitrary wave amplitude, the soliton dynamics 
in the presence of plasma inhomogeneity, the properties of linearly polarized solitons, 
and last but not least the stability of such structures (e. g., see discussions in Ref. 

imi)- 

In conclusion, we point out that the first macroscopic evidence of soliton formation 
in the interaction of an intense (10 19 W/cm 2 ) laser pulse with an underdense plasma 
has been reported in |28| . Long-lived, macroscopic bubble-like structures have been 
detected through the deflection that the associated electric charge separation causes 
in a proton probe beam These structures are interpreted as the remnants of a 
cloud of relativistic solitons generated in the plasma by the ultraintense laser pulse. 
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